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Abstract 



We consider an extended version of the pair tunneling model including interlayer single 
particle hopping (ISPH) as a complementary process to pair tunneling. The normal state 
gap, as found in cuprates, is taken to suppress the effective ISPH in conformity with the 
experimental observations, and this in turn enhances the pair tunneling process. The 
effective ISPH involves a probability factor P for which we consider two choices and 
provide phenomenological arguments in favour of them. We address the issue of bilayer 
splitting by calculating the spectral density function and corresponding photoemission 
intensity curves and show that our calculations conform with the absence of bilayer 
splitting observed in ARPES experiments on Bi2212. We have also studied the temper- 
ature variation of the superconducting gap and ratio of the superconducting gap to T c . 
Our results, obtained for both the choices of P, are reasonably in good agreement with 
those from experiments on cuprate superconductors. A linear T-dependent choice of P, 
however, yields a precise match to the experimantal data of the temperature varying 
superconducting gap. 
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I. Introduction 



In order to explain various unusual properties of high-T c superconductors (HTS) 
a number of phenomenological models have been proposed. One such model is the 
interlayer pair tunneling model (ILPT), originally proposed by P. W. Anderson 0. The 
model is based on the postulate that the coherent single particle tunneling between two 
Cu02 layers of HTS is blocked owing to strong electronic correlation effect, whereas 
tunneling of Cooper pairs between the layers is allowed. This derealization process of 
the pairs gives rise to a substantial enhancement of pairing in each layer and raises the 
transition temperature ||. An important outcome of the ILPT model is that it yields 
considerably higher T c for a bilayer system compared to that for a single layer system, 
even with a weak or moderate in-plane attractive interaction. Thus, the model can 
naturally explain higher values of T c in bi-layered or triple-layered systems compared to 
those of monolayer (in reality well separated layers) systems and obviates the need for 
a strong in-plane attractive interaction to achieve high T c . 

In the ILPT model of Anderson |], 0] the coherent single particle tunneling between 
layers is assumed to be completely blocked, even though the bare (interlayer) hopping 
rate, as obtained from the electronic structure calculation is substantial ||. The idea 
of absence of interlayer single particle tunneling originates from the anomalous c-axis 
transport properties of cuprates. For underdoped systems rapid increase in the c-axis 
resistivity with decreasing temperature and absence of a Drude peak in the optical 
conductivity along the c-axis || indicate absence of coherent single particle hopping 
between the layers 0. Anderson argued that these anomalies are due to the strong 
electronic correlations present in cuprates and are characteristics of a non-fermi liquid 
behaviour. However, for sufficiently overdoped cuprates the c— axis resistivity shows 
metallic like temperature dependence as that in the ab— plane and the Drude peak appars 
in the (c— axis) optical conductivity spectrum |5|, []]. These facts suggest presence of 
interlayer single particle hopping (ISPH) in case of overdoped systems. In fact, it is 
generally beleived that the confinement of the electrons in Cu02 planes is a characteristic 
feature of underdoped cuprates where correlation has its strongest effect. For overdoped 
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systems, on the other hand, the effect of correlation becomes weaker, hence the charge 
carriers are expected to behave more like a Fermi liquid. 

Recently, it has been established that a normal state pseudo gap (E g ) exists in high-T c 
systems |7], |j. This gap is momentum dependent and is given by 

E g (k) = -E g \Cosk x a — Cosk y a\ (1) 

where a is the in-plane lattice constant and the half factor is introduced so that the 
maximum value of the gap for any particular doping becomes E g . For underdoped 
systems E g is large, it decreases gradually with doping and becomes small or negligible 
for overdoped systems |7|, ||, §]. An interesting correlation has been observed between the 
normal state gap and the c-axis transport in cuprates. For underdoped system, where 
Eg is large, the interlayer single particle hopping is negligible, while, for overdoped 
system, where E g is small, it is allowed. Hence, a relationship, that the normal state 
gap is connected with the suppression of the interlayer single particle hopping, could be 
inferred. 

Motivated by the above observations and inference, recently we proposed an extended 



pair tunneling model |10| , |TT] involving both ISPH and interlayer pair tunneling [KJ, [11]], 
which is an extension of Anderson's ILPT model 0. We assume that the particles can 
tunnel between two layers via two channels: (i) single particle hopping and (ii) pair 
tunneling, and these two processes are complementary. We introduced a probability 
factor (P) for the ISPH which decreases very fast with the increase of E g [|Kj. As a 
result, for underdoped systems where E g is large, the ISPH is practically blocked and 
the only relevant interlayer coupling is through the pair tunneling. Our proposed model 
thus becomes equivalent to Anderson's ILPT model 0] for highly underdoped systems. 
For the overdoped case, where E g is small, the ISPH within our extended model becomes 
appreciable with the reduction of pair tunneling process. Thus, the extended model is 
consistent with the implications of the transport data which predict presence of ISPH 
for sufficiently overdoped systems. Certain important results of the extended ILPT 
model are: (i) it can naturally give high values of the ratio of the superconducting gap 
to T c , (ii) with increasing interlayer coupling T c may increase or decrease depending 



on the relative dominance of pair tunneling or ISPH, (iii) with the doping and other 
parameters remaining unchanged, T c becomes higher for materials with larger E g , (iv) 
temperature varying superconducting gap data is in closer agreement with that within 
the extended ILPT model than the original ILPT model. These results are similar to 
those characterizing the high-Tc cuprate superconductors. Some of these have already 
been presented elsewhere |TU|, |TT] . 



The presence of ISPH in our model may raise certain relevant questions regarding 
the issue of bilayer splitting in cuprate superconductors. The ISPH, in general, gives rise 
to splitting of bands which should be reflected in the electronic density of states (DOS) 
as well as in the spectral function, as the presence of a two peak struture. However, 
in recent angle resolved photoemission spectroscopy experiments (ARPES) JTB 



no evi- 



dence for bilayer splitting is observed in Bi2212 system. Therefore, it becomes important 
to investigate whether such observation is contradicted within our proposed extended 
model. In this communication we particularly focus on this issue. By a detailed study of 
the ARPES intensity curves, we show that the absence of bilayer splitting, as observed 
in experiments, can be understood within our proposed extended model. Our calcula- 



tions are carried out including realistic band structure of Bi2212 |L4| and the results are 
compared with the experimental data. 



In the proposed extended ILPT model |T0| , [TT| we consider the effective ISPH along 



with a probability factor P, which is suppressed by E g . Since the microscopic origin 
of Eg is still unknown and the proper functional dependence of P on E g is difficult to 
predict, we make two choices for P. In the first one |10| we consider 

P = e~ E ^ T (2) 



which is henceforth referred to as the exponential form. In another form [O 

T 



(3) 



E g + Te-^/ T 

which will be called as T -linear form. The exponential form could signify that the charge 
carriers has to overcome a gap for being available for the ISPH, whereas the T-linear 
form signifies no gap, but a simple power law dependence of P on T and E g . The latter 



form could find a justification within the spin-charge separation picture. Towards the 
end, we present arguments regarding these choices of P. 

In this communication we also make a comparative study of the properties related 
to the superconducting gap with the above mentioned two different choices of the prob- 
ability factor P. We find that, in general, the results are promising for both the choice, 
but the T-linear choice yields a better match to the experimental gap-variation data flTj) 
than the exponential choice. 



II. Formalism 



The model hamiltonian for the coupled bilayer system is given by 11 



where 



H = H N + H S (4) 



H N =Y, fa - ri c kMl + E {*f (*)c£M2 + h.c} (5) 



and 



Hs = - E c<3jc V$ + M 

- E {^(^cifc^a.cg + *.c} (6) 

The operator cjS + (cS) is the fermion creation (annihilation) operator with momen- 
tum k, spin | and in the layer i(= 1,2). if^v describes the band energy of the bilayer 
system, where e& is the band dispersion in a layer, t e f is the effective interlayer single 
particle hopping matrix element and \x is the chemical potential. Hs represents the in- 
teraction part of the hamiltonian and it leads to superconductivity in the system below 
the transition temperature T c . V^k' is the in-plane pairing interaction which is respon- 
sible for the formation of Cooper pairs in a layer. T£ s represents the effective interlayer 
pair tunneling matrix element. For the pair tunneling term we have considered only the 
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diagonal pair hopping , where quasi particle momentum is conserved during tunneling. 
Effective ISPH is taken to be tf(k) = t b ± (k)P, where t b ± (k) = t ± ((cosk x a - cosk y a)/2) 2 
is the fc-dependent ISPH as predicted by the band structure calculations U with t± 
being the bare ISPH matrix element, and P is the ISPH probability factor for which 
we consider the two forms as given in Eqs. (2) and (3). Pair tunneling, quite naturally, 
should be connected to the ISPH probability, since only those particles which are not 
taking part in ISPH, are available for the process. Tunneling of pairs, is a two particle 
process, hence, it occurs with a probability (1 — P) 2 |10|| . The effective pair tunneling 
matrix element is, thus, given by T^ s {k) = T p (k) (1 — P) 2 . The (bare) pair tunneling, 
following Ref. 0, is taken as T p (k) = [{t±{k)) 2 /\ti\], where t± is the nearest neighbour 
hopping matrix element of the a6-plane band dispersion. It may be noted that for large 
Eg or at very low temperatures the ISPH probability factor P becomes very small within 
our model and all the results predicted by our model become same as that of Anderson's 
ILPT model. In fact, in the limit E g /T — > oo (P — > 0), our model reduces exactly to 
the original ILPT model by Anderson. 

Mean field decoupling of four fermion terms in Hs yields 



Hs = -£ [A* + /».<:] (7) 

i,k 

where A k is the superconducting gap and is given by 

a, = A, fc = e <c V® ) + T f( fc ) < c -l4 } > ( g ) 

k' 

The layers i and j are equivalent, since by symmetry, the in-plane pairing average is 
identical in both the layers. 

It is well known that presence of interlayer single particle hopping between two layers 
produces bonding and antibonding bands. The corresponding annihilation operators are 
defined as cjj.j = (c^ — c^) and cffi = (c^ + cj^) and H N becomes diagonal in 
this new representation. The full hamiltonian may be written in terms of these new 
operators as 



H = E 



(9) 



where % is now the band index and % = — (+) represents the bonding (antibonding) 
band. Here, = (e k - //) ± tf(k) are the normal state band energies. From the 
hamiltonian of Eq.(9) the self consistent equations for the chemical potential and the 
superconducting gap are obtained as 



and 



k k 

E A ^M' (x(E k ,)+ X (E+))/2 

A fc = -r- (11) 

1- Tf(k) (x(E k ) + X (Et))/2 



where E k = \J + A|, are the quasiparticle energies in the superconducting state, 



x(E k ) = ^^tanh (^—^-J, ft = 1/T (in a scale of k B =l), 8 — 1 — n, with n being 
the number of electrons per site, and iV is the total number of lattice sites. The in- 
plane attractive pairing interaction is considered here only between nearest neighbors 
for which V^' is separable as = VrjkTjki, which also makes the /c-dependence of 
A/,, to be separable. In our calculations, the d x 2_ y 2 symmetry of the pairing state is 



considered because of growing evidence of the same in high-T c cuprates [16, 17]. This 
implies rj^. = (cos k x a — cosk y a)/2. 

Finally, the equation for the superconducting gap is obtained as 



1 i Vl {x{E k ) + X {Et)) 2 



AV N 



E r ( 12 ) 

* l-Tf{k)(x{E- k ) + x{Et))/2 

where the parameter V is the nearest neighbour in-plane attractive interaction strength. 

The ARPES intensity is proportional to the product of the spectral density function 
and the fermi distribution function. The spectral function in the superconducting state 
is given by 

A{k,uj) = A~(k,u) + A + (k 1 uj) (13) 

which includes contributions from both the (bonding and antibonding) quasiparticle 
bands. Spectral functions for the two bands are given by 



A ± (k,u) = - 



(uifv (vt) 2 r 



+ 



(u-Etf + T* {uj + Etf + T* 



(14) 
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i 2 4-2 i ( P 

where (u k ) = 1 — (v k ) = ~ ^ 1 + ) and T is a phenomenological linewidth parameter 



accounting for the broadening of the quasiparticle states due to finite lifetime. 

Considering the experimental energy resolution factor the intensity of the ARPES is 
given by [|i~8|l 

/oo 
R(uo - J)f(J)A(k',J)dJ (15) 

where R(uj) is the Gaussian energy resolution function and f(w) = \j{er* + 1) is the 
fermi distribution function. 



III. Results and Discussion 



A relevant point of contention, that may be raised regarding the validity of our model, 



as proposed in Ref. [10], [11], is that in experiments no evidence of bilayer splitting is 
observed in Bi2212. This suggests that the interlayer single particle hopping would be ab- 
sent, while a finite ISPH is present in our model. So, at first we present the results about 
the ARPES intensity within our model. By numerical solution of the self-consistent equa- 
tions (10) and (12) we evaluate the superconducting g £ip clS Qb function of temperature and 
hence determine the spectral density function as well as the ARPES intensity in the su- 
perconducting phase. For the in-plane band dispersion we consider a six parameter tight 
binding band structure, [to, t x , t 2 , h, t 4 , h] = [0.131,-0.149,0.041,-0.013,-0.014,0.013] 
eV, where to is the bare Wannier orbital energy, t\ nearest neighbor, t 2 next nearest 



neighbor etc. hopping matrix elements, as described |14j] and used in previous publica- 



tions |T9[. We take t± = 40meV and the in-plane attractive interaction V = 70meV, 
which are reasonable. To begin with, we consider the ARPES intensity and single par- 
ticle DOS for the exponential choice of P. 

The ARPES intensity I(k,oo) has been calculated using Eqs.(13)-(14) at the Brillouin 
zone point k = (ir, 0) for optimal doping and for different values of the energy resolution 
characterised by the full width at half maximum (FWHM) values of the Gaussian energy 
resolution function. The value of E g at optimal doping is taken from a best fit analysis 
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of the theoretical curve for the temperature variation of the superconducting gap to the 
experimental data, which is presented later. The plots of I(k,u) versus u at different 
temperatures are shown in Fig.l, where we take a very small value for the linewidth 
parameter T ~ lmeV. A perfect energy resolution is assumed in Fig. la. The figure 
shows that even within the limit of perfect resolution (FWHM=0) no bilayer splitting 
could be observed at T=13K, while at T=40K a splitting is observable in I(k,u>), but 
the difference in energy between the two peaks is only ~ 8meV . Fig. lb shows the 
evolution of I(k,u) at T=40K from a two peak structure to a single peak structure as 
the FWMH of the energy resolution function increases from 2meV to lOmeV. The 
experimental energy resolution for the ARPES, as quoted by different workers []i~8 , p0| 



is ~ 19meV, which is larger than that we have considered. Furthermore, at T = AO K 



the linewidth parameter in cuprates would be much larger than 1 meV [pi] , which will 
further broaden the spectral function and hence the I(k,u). So, it is clear that no 
bilayer splitting could be observed in the ARPES experiments within our model for the 
parameter space, relevant for high-T c cuprates, and for the exponential choice of P. 

Presence of the ISPH in a bilayer system would also produce a splitting in the normal 
state electronic DOS. So a study of the normal state DOS could provide a way to examine 
the significance of the bilayer splitting. In Fig. 2 we plot the electronic DOS at different 
temperatures and for different values of the linewidth parameters T. It may be noted 
that the effective ISPH in our model is temperature dependent and it makes the DOS 
to be dependent on temperature. Fig. 2a shows the plot of the DOS at T=10K. It is 
found that even with T=0, no evidence of bilayer splitting in the DOS could be observed 
within our model for the exponential choice of P at very low temperatures. For T= 
50K, bilayer splitting is quite evident for T=0. However, as the value of T is increased to 
14meV only one broad peak is observed in the DOS. For high-T c cuprates the value of 
T is quite large at high temperatures |2T| . Thus, our study shows that bilayer splitting 
may not be observable experimentally even if the normal state DOS is probed directly. 

In Fig. 3 we present the plot of the intensity of the ARPES (k = it, 0) for the T-linear 
form of P and for a value of the linewidth parameter T = lmeV. For this form of 
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P the suppression of the effective ISPH by E g at low temperatures is much less rapid 
than the exponential form of P. Consequently, at T = 13K the I(k,u) shows a two 
peak structure when the energy resolution is perfect (FWHM=0). However, even for 
FWHM as small as 3 meV the I(k,u) takes a single peak form at T = 13K. Fig. 3b 
shows how the two peak form of the I(k,u) at T = 40 K evolves into a single peak 
with decreasing energy resolution or increasing value of FWHM of the Gaussian energy 
resolution function. A single peak form for I(k, u) is obtained at T = AOK for a value of 
T as small as 1 meV and a value of FWHM of the energy resolution function as 12 meV, 



which is less than the FWHM values of the resolution function in experiments [18 



It should be noted that consideration of higher values of T would produce increased 
broadening in the shape of the I(k, u), hence a single peak form would be obtained even 
with a better energy-resolution (a lower value of FWHM) in actual experiments. 

In previous works |T(J [IT] we studied the superconducting phase diagram and some 
of the superconducting properties within our model for different E g values. The value 
of E g changes with doping in cuprates. However, for optimal doping E g has a fixed 
value for any particular cuprate material. In this communication we study different 
superconducting properties only at optimal doping and make a comparison of the results, 
obtained with two choices of P, with experimental data. In Fig.4 we present the best-fit 
curves for the temperature variation of the superconducting gap for two choices of P 
alongwith the experimental data points. The curve for the Anderson's pair tunneling 
model (Eg = oo) is also shown. It is clearly seen that predictions within our model 
are much closer to the experimental data than that with the original ILPT model. For 
the exponential choice of P, it is found that E g ~ 6 — 7meV gives a very good fit 
to the experimental data of Bi2212 at optimal doping. This value of E g is very close 
to the transition temperature of Bi2212 system at optimal doping. It is interesting to 
note that the value of E g ~ T c at optimal doping has been suggested in cuprates from 
experimental observations also 0. 

However, with the exponential form of P we find that it is difficult to obtain good fit 
to the experimental data points for the gap variation at low temperatures. A precise fit 
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to the experimental data is possible (Fig.4) for the second choice of P, where it increases 
almost linearly with temperature (Eq.(3)). However, the value of E g , required for such 
a fit is higher than that for the exponential choice of P. It should be noted that due to 
the absence of any microscopic derivation for the correct form of P, we make possible 
choices which are purely phenomenological. In the exponential choice, E g appears as a 
gap to hinder the charge transport along the c— axis. But, in the T-linear choice there 
is no direct gap and the conditions that P in Eq.(3) has to satisfy, are: (i) P increases 
linearly with T while decreases with E g for T <C E g and (ii) P = 1 for E g = 0. Clearly, 
these conditions could be satisfied even if a prefactor is added to E g . Thus the value of 
E g , as in Eq.(3), may differ from the actual value by a constant factor. 

In high-T c systems the ratio of the superconducting gap to T c is very high compared 
to the conventional superconductors. The value of the gap ratio (2A /T C ) in cuprates 



is reported to be within the range 6-9 | 5If . In Fig. 5 we plot the gap ratio as a function 
of interlayer single particle hopping matrix element (t±) for two different values of E g 
and for both the forms of P, for a doping level which is optimal for t± = AOmeV. The 
ratio increases with increasing value of t± . It is seen that both choices of P yield quite 
high values of the gap ratio, but for the T-linear choice of P the ratio is larger. For 
t± = AOmeV, which is realistic for layered cuprates, the value of the gap ratio lies within 
the range 5.0-7.5 (depending on E g and the choice of P). Thus, our model could yield 
high values for the gap ratio in agreement with the experimental findings in cuprates. 

It is observed that the gap ratio is higher for a smaller value of E g (Fig.5). Within 
our model a lower value of E g yields a lower value for T c m (transition temperature at 
optimal doping) [|1(]] and this is consistent with the experimental findings |9j . The results 
of Fig. 5, thus, predict a higher value of the gap ratio for a cuprate system with lower 
value of T c m . 

As found in the results presented above that the qualitative characteristics of them 
do not depend much on the specific form of the probability factor P used. While 
making a choice, one only needs to ensure that the probability factor P, and as a result 
the effective ISPH, follows the observed doping dependence of c-axis transport, that is 
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the effective c— axis hopping remains heavily suppressed in the underdoped side and at 
low temperatures. In experiments, E g is found to grow in magnitude |7|, || towards 
underdoing. Thus E g , on phenomenological grounds, could be taken to suppress c-axis 
hopping, and hence our choices of P could follow. 

It should be mentioned that, an exponential suppression of effective c-axis hopping 
towards underdoping was suggested || from the experimental results in Y-123 system 
where a metallic like behavior of c-axis resistivity at high temperatures and a semicon- 
ducting behavior at low temperatures have been observed. But, our calculations show 
that the T-linear choice of P gives a better match to the experimental gap-variation 
data. One could probably search for a realization of this T-linear choice within within 
the RVB picture ]E| . According to the RVB model, holons and spinons are the quasipar- 
ticles in a layer of strongly correlated system and a holon has to combine with a spinon 
to form a hole which can then hop from one layer to another. Consequently, c— axis 
hopping is proportional to the spinon density, which increases linearly with T within 
the RVB theory and hence, the effective c— axis single particle hopping is expected to 
be proportional to T. 



IV. Conclusions 



Investigation of the spectral density function and corresponding intensity curves for the 
ARPES within our proposed extendeded pair tunneling model shows that the model 
has no contradiction with the absence of bilayer splitting observed in Bi2212 by ARPES 



experiment |L3J . Studies with two different forms for the probability factor P, associated 
with the effective interlayer single particle hopping, show that both the forms could 
yield high values for the ratio of the superconducting gap to T c , as observed in cuprates. 
Temperature varying superconducting gap from our calculations with both the forms 
agrees well with the experimental data and the agreement is much better compared to 
that within the original ILPT model. A precise fit to the experimental gap variation 
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data is achieved with the T-linear choice of P. Another prediction of our model is that 
a system which has a lower value of E g should have a lower value of T c m and a larger 
value of the gap ratio. 

It may be mentioned that Anderson's ILPT model is applicable in the underdoped 
region where non-fermi liquid behavior dominates and the ISPH is negligible, whereas 
our extended model has the flexibility of being applicable in a wide region from the 
underdoped to a highly overdoped regime where the ISPH is appreciable. 

In the model, presented here, we have considered the effect of E g only on the interlayer 
hopping. However, it seems natural to think that, E g should also affect the in-plane 
charge dynamics. In this connection, it has been suggested that E g reduces the in-plane 
electronic density of states || . It would be interesting to study the extended model with 
the pseudogap E g affecting both the in-plane and out of plane charge dynamics. This 
would constitute future communications. 
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Figure captions: 

Fig.l. Plot of the ARPES intensity I(k,u>) in arbitrary units versus uo (eV) for the 
exponential choice of P. Different values of FWHM of the Gaussian energy resolution 
function and temperatures are as given in the figure. The linewidth parameter is taken 
to be T = 1 meV. 

Fig. 2. Plot of the electronic density of states iV(£) versus £ for the exponential choice 
of P. Temperatures and different values of linewidth parameter T are mentioned in the 
figure. 

Fig. 3. The ARPES intensity I(k,u) versus u for the T-linear choice of P with the 
linewidth parameter T = 1 meV. Values of FWHM and temperatures are listed in the 
figure. 

Fig. 4. Finite temperature gap scaled to its zero temperature value (A™ ax (T) / A™ ax (0)) , 
as a function of reduced temperature (T/T c ). Solid line is that in the Anderson limit 
(Eg = oo), solid square symbols are experimental data from Ref.[[0J and dashed lines 
are from our calculations for two different forms of the probability factor P. 

Fig. 5. Maximum value of the gap-ratio (2A™ a:r (0)/T c ) versus bare interlayer coupling 
(t±), for different E g values and different choices of P, as mentioned in the figure. 
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